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Abstract 

When independent Bose-Einstein condensates (BEC), described quantum mechanically by Fock 
(number) states, are sent into interferometers, the measurement of the output port at which the 
particles are detected provides a binary measurement, with two possible results ±1. With two inter- 
ferometers and two BEC's, the parity (product of all results obtained at each interferometer) has all 
the features of an Einstein-Podolsky-Rosen quantity, with perfect correlations predicted by quantum 
mechanics when the settings (phase shifts of the interferometers) are the same. When they are differ- 
ent, significant violations of Bell inequalities are obtained. These violations do not tend to zero when 
the number A'' of particles increases, and can therefore be obtained with arbitrarily large systems, but 
a condition is that all particles should be detected. We discuss the general experimental requirements 
for observing such effects, the necessary detection of all particles in correlation, the role of the pixels 
of the CCD detectors, and that of the alignments of the interferometers in terms of matching of the 
wave fronts of the sources in the detection regions. 

Another scheme involving three interferometers and three BEC's is discussed; it leads to Greenberger- 
Horne-Zeilinger (GHZ) sign contradictions, as in the usual GHZ case with three particles, but for an 
arbitrarily large number of them. Finally, generalizations of the Hardy impossibilities to an arbitrarily 
large number of particles are introduced. BEC's provide a large versality for observing violations of 
local realism in a variety of experimental arrangements. 

PACS numbers: 03.65.Ud, 03.75.Gg, 42.50.Xa 

The original Einstein-Podolsky-Rosen (EPR) argument [1] considers a system of two microscopic 
particles that are correlated; assuming that various types of measurements are performed on this system 
in remote locations, and using local realism, it shows that the system contains more "elements of reality" 
than those contained in quantum mechanics. Bohr gave a refutation of the argument [2] by pointing 
out that intrinsic physical properties should not be attributed to microscopic systems, independently of 
their measurement apparatuses; in his view of quantum mechanics (often called "orthodox"), the notion of 
reality introduced by EPR is inappropriate. Later, Bell extended the EPR argument and used inequalities 
to show that local realism and quantum mechanics may sometimes lead to contradictory predictions [3] . 
Using pairs of correlated photons emitted in a cascade, Clauser et al. [5] checked that, even in this case, 
the results of quantum mechanics are correct; other experiments leading to the same conclusion were 
performed by Fry et al. [6], Aspect et al. [7], and many others. The body of all results is now such that 
it is generally agreed that violations of the Bell inequalities do occur in Nature, even if experiments are 
never perfect and if "loopholes" (such as sample bias [8-10]) can still be invoked in principle. All these 
experiments were made with a small number of particles, generally a pair of photons, so that Bohr's point 
of view directly applies to them. 



1 



In this article, as in [11] wc consider systems made of an arbitrarily large number of particles, and 
study some of their variables that can lead to an EPR argument and Bell inequalities. Mermin [12] has 
also considered a physical system made of many particles with spins, assuming that the initial state is a 
so called GHZ state [13,14]; another many-particle quantum state has been studied by Drummond [15]. 
Nevertheless, it turns out that considering a double Fock state (DFS) with spins, instead of these states, 
sheds interesting new light on the Einstein-Bohr debate. The reason is that, in this case, the EPR elements 
of reality can be macroscopic, for instance the total angular momentum (or magnetization) contained in 
a large region of space; even if not measured, such macroscopic quantities presumably possess physical 
reality, which gives even more strength to the EPR argument. Moreover, one can no longer invoke the 
huge difference of scales between the measured properties and the measurement apparatuses, and Bohr's 
refutation becomes less plausible. 

Double Fock states with spins also lead to violations of the Bell inequalities [16,17], so that they 
are appropriate for experimental tests of quantum violations of local realism. A difficulty, nevertheless, 
is that the violations require that all N spins be measured in N different regions of space, which may 
be very difficult experimentally if N exceeds 2 or 3; with present experimental techniques, the schemes 
discussed in [16, 17] are therefore probably more thought experiments than realistic possibilities. Here we 
come closer to experiments by studying schemes involving only individual position measurement of the 
particles, without any necessity of accurate localization. 

With Bose condensed gases of metastable helium atoms, micro-channel plates indeed allow one to 
detect atoms one by one [18, 19]. The first idea that then comes to mind is to consider the interference 
pattern created by a DFS, a situation that has been investigated theoretically by several authors [20-24], 
and observed experimentally [25]. The quantum effects occurring in the detection of the fringes have 
been studied in [26,27], in particular the quantum fluctuations of the fringe amplitude; see also [28] for a 
discussion of fringes observed with three condensates, in correlation with populations oscillations. But, 
for obtaining quantum violations of Bell type inequalities, continuous position measurements are not 
necessarily optimal; it is more natural to consider measurement apparatuses with a dichotomic result, 
such as interferometers with two outputs, as in [29,30]. Experimentally, laser atomic fluorescence may 
be used to determine at which output of an interferometer atoms are found, without requiring a very 
accurate localization of their position; in fact, since this measurement process has a small effect on the 
measured quantity (the position of the atom in one of the arms), one obtains in this way a quantum 
non-demolition scheme where many fluorescence cycles can be used to ensure good efficiency. 

Quantum effects taking place in measurements performed with interferometers with 2 input and 2 
output ports (Mach-Zhender interferometers) have been studied by several authors; refs [31-33] discuss the 
effect of quantum noise on an accurate measurement of phase, and compare the feeding of interferometers 
with various quantum states; refs [34-36] give a detailed treatment of the Heisenberg limit as well as of 
the role of the Fisher information and of the Cramer-Rao lower bound in this problem. But, to our 
knowledge, none of these studies leads to violations of Bell inequalities and local realism. Here, we will 
consider interferometers with 4 input ports and 4 output ports, in which a DFS is used to feed two of 
the inputs (the others receive vacuum), and 4 detectors count the particles at the 4 output ports - see 
Fig. [Tl we will also consider a similar 6 input-6 output case. This can be seen as a generalization of 
the work described by Yurke and Stoler in refs. [37,38], and also to some extent of the Rarity- Tapster 
experiment [39,40] (even if, in that case, the two photons were not emitted by independent sources). 

Another aspect of the present work is to address the question raised long ago by Anderson [41] in the 
context of a thought experiment and, more recently, by Leggett and Sols [42,43]: "Do supcrfluids that 
have never seen each other have a well-defined relative phase?" . A positive answer occurs in the usual 
view: when spontaneous symmetry breaking takes place at the Bose- Einstein transition, each condensate 
acquires a well-defined phase, though with a completely random value. However, in quantum mechanics, 
the Bose-Einstein condensates of superfluids are naturally described by Fock states, for which the phase 
of the system is completely undetermined, in contradiction with this view. Nevertheless, the authors of 
refs. [20-23] and [29,44] have shown how repeated quantum measurements of the relative phase of two 
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Fock states make a well-defined value emerge spontaneously, with a random value. This seems to remove 
the contradiction; considering that the relative phase appears under the effect of spontaneous symmetry 
breaking, as soon as the EEC's are formed, or later, under the effect of measurements, then appears as 
a matter of personal preference. 

But a closer examination of the problem shows that this is not always true [16,27,42,43]: situations do 
exist where the two points of view are not equivalent, and where the predictions of quantum mechanics 
for an ensemble of measurements are at variance with those obtained from a classical average over a 
phase. This is not so surprising after all: the idea of a pre-existing phase is very similar to the notion of 
an EPR "element of reality" - for a double Fock state, the relative phase is nothing but what is often 
called a "hidden variable". The tools offered by the Bell theorem are therefore appropriate to exhibit 
contradictions between the notion of a pre-existing phase and the predictions of quantum mechanics. 
Indeed, we will obtain violations of the BCHSH inequalities [45], new GHZ contradictions [13,14] as well 
as Hardy impossibilities [46,47]. Fock-state condensates appear as remarkably versatile, able to create 
violations that usually require elaborate entangled wave functions, and produce new iV-body violations. 

A preliminary short version of this work has been published in [48] . The present article gives more 
details and focusses on some issues that will inevitably appear in the planning of an experiment, such 
as the effect of non-perfect detection efficiency, losses, or the geometry of the wavefronts in the region 
of the detectors. In § [T] we basically use the same method as in [48] (unitary transformations of creation 
operators), following refs [37] and [38], but include the treatment of losses; in § [51 we develop a more 
elaborate theory of many-particle detection and high order correlation signals, performing a calculation in 
3A'^ configuration space, and including a treatment of the geometrical effects of wavefronts in the detection 
regions (this section may be skipped by the reader who is not interested in experimental considerations); 
finally, §[3] applies these calculations to three situations: BCHSH inequality violations with two sources, 
GHZ contradictions with three sources, and Hardy contradictions. Appendix I summarizes some useful 
technical calculations; appendix II extends the calculations to initial states other than the double Fock 
state ((!]), in particular coherent and phase states. 

1 Quantum calculation 

We first calculate the prediction of quantum mechanics for the experiment that is shown schematically in 
Fig. dJ Each of two Bose-Einstein condensates, described by Fock states with populations Na and Nj^, 
crosses a beam splitter; both are then made to interfere at two other beam splitters, sitting in remote 
regions of space Da and Db- There, two operators, Alice and Bob, count the number of particles that 
emerge from outputs 1 and 2 for Alice, outputs 3 and 4 for Bob. By convention, channels 1 and 3 are 
ascribed a result ?7 = +1, channels 2 and 4 a result 77 = —1. We call rrij the number of particles that 
are detected at output j (with j = 1, 2, 3, 4), uia = rni + m2 the total number of particles detected by 
Alice, niB = + 7714 the number of particles detected by Bob, and M = niA + the total number 
of detected particles. From the series of results that they obtain in each run of the experiment, both 
operators can calculate various functions A{r]i, ■■rjmA) and B{r/„iA+ij ■■Vm) of their results; we will focus 
on the case where they choose the parity, given by the product of all their 77's: A — (—1)™^ for Alice, 
B = (—1)"'* for Bob; for a discussion of other possible choices, see [17]. 

We now calculate the probability of any sequence of results with the same approach as in [48] . This 
provides correct results if one assumes that the experiment is perfect; a more elaborate approach is 
necessary to study the effects of experimental imperfections, and will be given in § [51 
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Figure 1: Two Fock states, with populations Na and N^, enter beam splitters, and are then made to 

interfere in two different regions of space Da and Db, with detectors 1 and 2 in the former, 3 and 4 in 
the latter. The number of particles nij in each of the channels j = 1, 2, 3, 4 are counted. 



1.1 Probabilities of the various results 

We consider spinless particles and we assume that the initial state is: 



\^o) = \Na,Np) = 



K)^ |o) 



(1) 



where |0) is the vacuum state; single particle state a corresponds to that populated by the first source, 
(} to that populated by the second source. The destruction operators ai • • • 04 of the output modes can 
be written in terms of those of the modes at the sources Oq, a^, and aai (including the vacimm input 
modes, aa' and a^', which are included to maintain unitarity) by tracing back from the detectors to 
the sources, providing a phase shift of it/ 2 at each reflection and C or ^ at the shifters, and a l/\/2 for 
normalization at each beam splitter. Thus we find: 
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Since a^' and a^' do not contribute we can write simply: 



ai 

as 
04 



— e"»aa + ap 



In short, we write these expressions as: 



(3) 



(4) 



We suppose that Alice finds mi positive results and m2 negative results for a total of m,A measurements; 
Bob finds ma positive and m4 negative results in his itib total measurements. The quantum probability 
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of this scries of results is the squared modulus of the scalar product of state |<i>o) by the state associated 
with the measurement: 

V{m^,m2,7m,m^) = , ] , , | (0| (ai)"^ • • • (04^ |iV„, TV^^) |' (5) 
mi]m2'-m3\m4\ 

where the matrix element is non-zero only if: 

mi + 1712 + my, + 1714 = Na + N(3 = N (6) 

We can calculate this matrix element by substituting ([1]) and (j4]) and expanding in binomial series: 
(0| (ai)'- • • • (ai)™^ \N^,Np) = \ (0| f[ {v,^a^ + v^pap^ (4)'"° U) |0) 

where pi3j — rrij — paj for any j. The matrix element at the end of this expression is: 

But by definition the sum of all p's is equal to mi + m2 + TO3 + TO4 which, according to ([H]), is A^q, + Np; 
the two Kronecker delta's in Q are therefore redundant. For the matrix element to be non zero and equal 
to NalNpl, it is sufficient that the difference between the sums pai + • • • + PaA and ppi + • • • + pp^ be 
equal to Na — Nf), a condition which we can express through the integral: 

^ i{N0~N^+p^l + ---+ p^4~p01 + --- -^P,34)M - X.. „ ^ ^ (Q) 

~ '^Na-Nff, PoiH h pai-ppi —P/34, 

When this is inserted into ([7]), in the second line every Vja becomes Vjae^^ , every Vjp becomes Vjf^e~^^, 
and we can redo the sums and write the probability amplitude as: 

x/raV r ge'^"''"""^'' n (-."^''" + ^^P^-'T (10) 
Thus the probability is: 

V {mi, 7712, 1713, rrii) ^ j j j 7 ^ ^ II ^^j M ^ (W (H) 

mi\m2\m3\m,i\ J 2tt J ^ 2tt j-J-^^ 

with: 

nJ{^^) =vjae"' + Vjpe-'>' (12) 
Each of the factors can now be simplified according to: 

which, when Q is inserted, gives: 

- [cos (/i — ^') ± COS (C + /-i + /^')] O'' 2 ^P ^ ^'^ ^'^^ + /-i + /u')] 
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depending on the value of j. Now, if we define: 

A = /i + /i' 
A = fi- fi' 

we finally obtain: 

\m2lm3\m4\ 2tt 27r 

X [cos A + cos (C + A)]™^ [cos A - cos (C + A)]"^ [cos A + cos {6 - A)]"' [cos A - cos {6 - A)] 

where we have used A parity to reduce cosine. 



(14) 



mi\m2\ms\mi\ J_^2ttJ_^2tt (15) 



1.2 Effects of particle losses 

We now study cases where losses occur in the experiment; some particles emitted by the sources are 
missed by the detectors sitting at the four output ports. The total number of particles they detect is M, 
with M < A^; an analogous situation was already considered in the context of spin measurements [16,17]. 
We first focus on losses taking place near the sources of particles, then on those in the detection regions. 

1.2.1 Losses at the sources 

As a first simple model for treating losses, we consider the experimental configuration shown in Fig. [21 
where additional beam splitters divert some particles before they reach the input of the interferometer. 
If T and R arc the transmission and the refiection coefficients of the additional beam splitters, with: 



the unitary transformations become: 



R + T^l (16) 



«3 = ^ K + (17) 

05 = i\/R [aa] 
ae = i^/R [ap] 



(18) 



The probability amplitude associated with a series of results mi, ...ms, mg is now: 

(0| (ai)"^ ■ ■ ■ ja^r' ja^r' K)"^ Wa, Np) 
\/mi\....m5\ me! 

or, taking into account the last two equations ([T7]): 



/ Nal ^ Np\ (0| (ai)'"^ ■ ■ ■ (04)"" \Na - m^^Np - mg) 

\l {Na - m^y. m^l ^ [Np - me)l mel 

The fraction on the right of this expression can be obtained from the calculations of § 11.11 by just replacing 
Na by Na — ms, Np by Np — mg in (fTO| . With this substitution, the numerical factor in front of that 
expression combines with that of to give a prefactor: 



i?'"^-"--,/^ (20) 
V mslmg! 
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Figure 2: The experiment is the same as in figure [TJ but now we assume that two beam sphtters are 
inserted between the two sources and the inputs of the interferometer. Then, the total number of particles 
M measured at the output of the interferometer may be less than N . 



The next step is to sum the probabilities over and mg, keeping the four mi, ..m4 constant; we 
vary ms and mg with a constant sum N — M, where M is defined as: 

M = mi + m2 + ms + m4 < iV (21) 

The cos[{Ni3 — Na) A] inside the integral of p5|) arose from an exponential e*t'^''~^°^'^l , which now 
becomes: 

^i{Ni3-Na+m:i-m6)A ^22) 

SO that the summation over and mg reconstructs a power of a binomial: 

When the powers of R and T are included as well as the factors 2^^^ and 2^^"^^), equation (fTSjl is now 
replaced by: 

■mi\m2\m'i\m4 {N - M)\ J _^ 2t: J _^ 2t: 

X [cos A + cos (C + A)]™^ [cos A - cos (C + A)]"^ 
X [cos A + cos {e - A)]'"^' [cos A - cos {6 - A)]"* 

(24) 

where M is defined in (pTjl . This result is similar to ([T5)l , but includes a power of cos A inside the integral, 
which we will discuss in § 11.31 We note that this power of cos A introduces exactly the same factor as 
that already obtained in [16], in the context of spin condensates and particles missed in transverse spin 
measurements . 

If T = 1 and i? = 0, expression (|24|l vanishes unless M has its maximal value M = N\ then expression 
is recovered, as expected. If R and T have intermediate values, M has a probability distribution 
including any value less than iV, with of course smaller values favored when T is small and R large. 



VM{mi,m2,m3,m4) ^ 1 -2" —t^j 777^/ TTI cos [(iV^ - iVQ)A] [cosA] 



1.2.2 Losses at the detectors 

Instead of inserting additional beam splitters just after the sources, we can put them just before the 
detectors, as in Fig. [31 this provides a model for losses corresponding to imperfect detectors with quantum 
efficiencies less than 100%. 
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Figure 3: The experiment is the same as that in figure 2, but now 4 beam sphtters are inserted just 
before the 4 particle detectors, which sit in channels 1,2,3,4; the other channels, 1', 2', 3', 4' contain no 
detector. This provides a model for calculating the effect of limited quantum efficiencies of the detectors. 



Instead of 6 destruction operators as in ([T7)) . we now have 8; each operator aj {j = 1,2,3,4) 
corresponding to one of the 4 the detectors is now associated with a second operator a'^ corresponding to 
the other output port with no detector. For instance, for j ~ I, one has: 



ai = [e"^aa + ; a[ = — ^ [e*'^a„ + ap] (25) 

and similar results for j ^ 2,3,4. The calculation of the probability associated with results mi, ..7714 (with 
their sum equal to M) and m[, ..m'^ (with their sum equal to — AI) is very similar to that of § 11.11 
formula ([5]) becomes: 

V{m,, ..m,; m[, ..m',) = , , / , (0| (ai)"^ • • • (04)"^ x {aX'^ ■ ■ ■ {a^'' Wc^.N^) ' 

mi!..!m4! x m\\..\m^\ 

Since aj and a'j are almost the same operator (they just differ by a coefficient), the result is still given 
by the right hand side of (fT5|) . with the following changes: 

(i) each mj is now replaced by the sum nij + m'j 

(ii) mi!..!m4! in the denominator is multiplied by mi!..!m4! 
(ii) a factor x R^^^^ appears in front of the expression. 

Now, we consider the observed results mi, ..7114 as fixed, and add the probabilities associated with all 
possible non-observed values m'l, ..m4; this amounts to distributing N — M unobserved particles in any 
possible way among all output channels without detectors. The summation is made in two steps: 

(i) summations over m'^ and m'j at constant sum m[ + m'j = m^, and over mg and m'^ at constant 
sum m'^ + m'4 = m'g 

(ii) summation over m^ and m'g at constant sum m^ + m'g — N — M 
The first summation provides: 



,1 [cos A + cos (C + A)]™'i [cos A - cos (C + A)]"^ = [2cos A]™^ (26) 



E 
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and similarly for the summation over m'^ and m^. Then the last summation provides: 



y ,] , , [2 cos A]"^ [2 cos A]"'« = [4 cos A]"^-"^' (27) 

m'j^+m'^ = N-M ^ ^ 

At the end, as in § 11.2.11 we see that each unobserved particle introduces a factor cos A so that, in the 
integral giving the probability, a factor [cosA]''^"*^ appears; actually, we end up with an expression that 
is again exactly ((24|) . The presence of this factor inside the integral seems to be a robust property of 
the effects of imperfect measurements (for brevity, we do not prove the generality of this statement, for 
instance by studying the effect of additional beam splitters that are inserted elsewhere, for instance in 
other parts of the interferometer). 

1.3 Discussion 

The discussion of the physical content of equation (|15p and is somewhat similar to that for spin 
measurements [16]. One difference is that, here, we consider that the total number of measurements 
TUA = mi + 7712 made by Alice, as well as the total number of measurements ttib — 7713 + 7774 made 
by Bob, are left to fluctuate with a constant sum M; the particles emitted by the sources may localize 
in any of the four detection regions. With spins, the numbers of detections depends on the number of 
spin apparatuses used by Alice and Bob, so that it was more natural to assume that ttia and rriB are 
fixed. This changes the normalization and the probabilities, but not the dependence on the experimental 
parameters, which is given by ()15|) and (I24p . A discussion of the normalization integrals is given in the 
Appendix. 

1.3.1 Effects of Na, Np and of the number of measurements 

If the numbers of particles in the sources are different {Na 7^ ^p), a term in cos [(A^/j — A"q)A] appears 
in both equations; let us assume for simplicity that all N particles are measured (M = A^), so that 
equation (fTS]) applies, and for instance that A^q, > A''^. Then, in the product of factors inside the integral, 
only some terms can provide a non-zero contribution after integration over A; we must choose at least 
Na — Np factors contributing through cos (A), and thus at most N—{Na — Np) = 2Np factors contributing 
through the 9 dependent terms. Therefore 2A^^ is the maximum number of particles providing results 
that depends on the settings of the interferometer; all the others have equal probabilities 1/2, whatever 
the phase shift is. This is physically understandable, since (A'q, — A^^) particles from the first source 
unmatched particles from the other, and can thus not contribute to an interference effect. All particles 
can contribute coherently to the interference only if A^^ — Np. 

If the numbers of particles in the sources are equal {N^ = Np), the sources are optimal; equation 
contains the effect of missing some particles in the measurements. If the number of experiments M is 
much less than a very large A^, because cos A^^*^ peaks up sharpl}0 at A = 0, the result simplifies into: 

VM{mi,m2,m3,m,) , , , , / ^ [1 + cos (C + A)]™^! " cos (C + A)]™^ 

X [1 + cos {9 ~ A)]'"-^ [1 - cos {9 - A)]™* (28) 

We then recover "classical" results, similar to those of refs. [29] or [49]. Suppose that we introduce a 
classical phase A and calculate classically the interference effects at both beam splitters. This leads to 
intensities proportional to [1 + cos + A)] and [1 — cos {( + A)] on both sides of the interferometer in 

^Herc wc take the point of view where the A integration domain is between — 7r/2 and +7r/2; otherwise, we should also 
take into account a peak around A = tt. 
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Da, and similar results for Db- Now, if we assume that each particle reaching the beam splitter has 
crossing and reflecting probabilities that are proportional to these intensities (we treat each of these 
individual processes as independent), and if we consider that this classical phase is completely unknown, 
an average over 2n then reconstructs exactly (1^ . In this case, the classical image of a pre-existing phase 
leads to predictions that are the same as those of quantum mechanics; this phase will take a completely 
random value for each realization of the experiment, with for instance no way to force it to take related 
values in two successive runs. All this fits well within the concept of the Anderson phase, originating from 
spontaneous symmetry breaking at the phase transition (Bose-Einstein condensation): at this transition 
point, the quantum system "chooses" a phase, which takes a completely random value, and then plays 
the role of a classical variable (in the limit of very large systems). 

On the other hand, if TV — M vanishes, the peaking effect of cos A^~*^ does not occur anymore, A can 
take values close to 7r/2 , so that the terms in the product inside the integral are no longer necessarily 
positive; an interpretation in terms of classical probabilities then becomes impossible. In these cases, the 
phase does not behave as a semi-classical variable, but retains a strong quantum character; the variable 
A controls the amount of quantum effects. It is therefore natural to call A the "quantum angle" and A 
the "classical phase" . 

One could object that, if expression contains negative factors, this does not prove that the same 
probabilities V can not be obtained with another mathematical expression without negative probabilities. 
To show that this is indeed impossible, we have to resort to a more general theorem, the Bell/BCHSH 
theorem, which proves it in a completely general way; this is what we do in § [3] 



1.3.2 Perfect correlations 



We now show that, if Na and iV^ are equal and if the number of measurements is maximal (Af = N), 
when Alice and Bob choose opposite^ phase shifts (6 — they always measure the same parity. In this 
case, the integrand of (fTK|) becomes: 



[cos A -I- cos (A -t-C)] 



mi+m3 



COS A — COS (A -I- ()] 



m2+m4 



(29) 



which can also be written as: 



sin(A' + ^) 



sin(A" - ^) 



cos(A 



-1 m2+m4 



cos(A" - |) 



m2+™4 



with the following change of integration variable^: 



A' = A±A 



A" = A^ 



(30) 



(31) 



If, for instance, nii+m^ is odd, instead of A' one can take (A'-|-^/2) as an integration variable, and one can 
see that the integral vanishes because its periodicity - the same is true of course for the A" integration, 
which also vanishes. Similarly, if m2 -\- nii is odd, one can take (A' — and (A" — ^) as integration 
variables, and the result vanishes again. Finally, the probability is non-zero only if both mi -I- m^, and 
7712 + are even; the conclusion is that Alice and Bob always observe the same parity for their results. 
This perfect correlation is useful for applying the EPR reasoning to parities. 



^Usually, perfect correlations are obtained wlien tlie two settings are the same, not opposite. But, with the geometry 
shown in figure^ introduces a phase delay of source a with respect to source /3, while 9 does the opposite by delaying 
source /3 with respect to source a. Therefore, the dephasing effects of the two delays are the same in both regions Da and 
Dg when 9 = — f . 

^Using the periodicity of the integrand, one can give to both integration variables A' and A" a range [— tt, +7r]; this 
doubles the integration domain, but this doubling is cancelled by a factor 1/2 introduced by the Jacobian. 
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2 A more elaborate calculation 



The advantage of measuring the positions of particles after a beam sphtter, that is interference effects 
providing dichotomic resuhs, is that one has a device that is close to quantum non-demolition experiments. 
With a resonant laser, one can make an atom fluoresce and emit many photons, without transferring the 
atom from one arm of the interferometer to the other. This is clearly important in experiments where, as 
we have seen, all the atoms must be detected to obtain quantum non-local effects. On the other hand, it 
is well know experimentally that a difficulty with interferometry is the alignment of the devices in order 
to obtain an almost perfect matching of the wave front. We discuss this problem now. 

A general assumption behind the calculation of §[T]is that, in each region of space (for instance at the 
inputs, or at the 4 outputs), only one mode of the field is populated (only one operator is introduced 
per region). The advantage of this approach is its simplicity, but it nevertheless eludes some interesting 
questions. For instance, suppose that the energies of the particles emerging from each source differ by 
some arbitrarily small quantity; after crossing all beam splitters, they would reach the detection regions 
in two orthogonal modes, so that the probability of presence would be the sum of the corresponding 
probabilities, without any interference term. On the other hand, all interesting effects obtained in [48] 
are precisely interference effects arising because, in each detection region, there is no way to tell which 
source emitted the detected particles. Does this mean that these effects disappear as soon as the sources 
are not strictly identical, so that the quantum interferences will never be observable in practice? 

To answer this kind of question, here we will develop a more detailed theory of the detection of many 
particles in coincidence, somewhat similar to Glauber's theory of photon coincidences [50]; nevertheless, 
while in that theory only the initial value of the n-th time derivative was calculated, here we study the 
whole time dependence of the correlation function. Our result will be that, provided the interferometers 
and detectors are properly aligned, it is the detection process that restores the interesting quantum effects, 
even if the sources do not emit perfectly identical wavefronts in the detection regions, as was assumed 
in the calculations of § [TJ As a consequence, the reader who is not interested in experimental limitations 
may skip this section and proceed directly to § [31 

In this section we change the definition of the single particle states: a now corresponds to a state for 
which the wave function originates from the first source, is split into two beams when reaching the first 
beam splitter, and into two beams again when it reaches the beam splitters associated with the regions 
of measurement Da and Db', the same is true for state /?. In this point of view, all the propagation in 
the interferometers is already included in the states. We note in passing that the evolution associated 
with the beam splitters is unitary; the states a and P therefore remain perfectly orthogonal, even if they 
overlap in some regions of space. Having changed the definition of the single particle states, we keep ([T]) 
to define the N particle state of the system. 

2.1 Pixels as independent detectors 

We model the detectors sitting after the beam splitters by assuming that they are the juxtaposition of a 
large number Q of independent pixels, which we treat as independent detectors. This does not mean that 
the positions of the impact of all particles are necessarily registered in the experiment; our calculations 
still apply if, for instance, only the total number of impact in each channel is recorded. The only thing 
we assume is that the detection of particles in different points leads to orthogonal states of some part 
of the apparatus (or the environment), so that we can add the probabilities of the events corresponding 
to different orthogonal states; whether or not the information differentiating these states is recorded in 
practice does not matter. 

If the number of pixels Q is very large, the probability of detection of two bosons at the same pixel 
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is negligible. If we note N — Na + Np, this probability is boundc43 by: 

1 2 N N{N+1) 



(32) 



Q Q Q 2Q 

so that we will assume that: 

Q > (33) 

Moreover, we consider events where all N particles are detected (the probabilities of events where some 
particles are missed can be obtained from the probabilities of these events in a second step, as in [16]). 

2.2 Flux of probability at the pixels in a stationary state 

Each pixel j is considered as defining a region of space Aj in which the particles are converted into a 
macroscopic electric current, as in a photoniultiplier. The particles enter this region through the front 
surface Sj of the pixel; all particles crossing Sj disappear in a conversion process that is assumed to have 
100% efficiency. What we need, then, is to calculate the flux of particles entering the front surfaces of the 
pixels. It is convenient to reason in the 3iV dimension configuration space, in which the hyper-volume 
associated with the N different pixels is: 

Va, = Ai A2 A3.... Ajv (34) 

which has an front surface given by: 

5jv = 51 «) A2 (8) A3.... «) Aat + Ai 5*2 «) A3.... (g) Aat + .. + Ai (g) A2 «) A3.... (g) S'jv (35) 

The density of probability in this space is defined in terms of the boson field operator ^'(r) as: 

pjv(ri,r2,...,rA,) =^'t(ri)*(ri) *t(i.2)^-(r2) ... {r n)'^' {r n) (36) 

where we assume that all r^'s are different (all pixels are disjoint). The components of the 3A^ dimension 
current operator 3n are: 

f 2^ [*^(ri)V^'(ri)- V*t(ri)4-(ri)] *t(i.2)5.(r2) ...*t(r^)^(rjv)+ 

+ 2^*^(ri)*(ri) [*^(r2)V^'(r2)- V^'t(r2)4'(r2)] ...«'t(r^)*(r^)+ ^^^^ 

+ 2ls*^(i-i)*(i'i) *^(r2)*(r2) ... [^'t(rw)V*(r^)- V*t(r^)^(rjv)] 
In the Heisenberg picture, the quantum operator pjv obeys the evolution equation: 

^PAr(ri,r2, ...,rAr;i) + Vat • Jat = (38) 

where V^v is the N dimensional divergence. The flux of probability entering the 3A^ dimension volume 
Sn it then: 

.F(Ai, A2, ...An) = <^o\F (Ai) x G (A2) x ....G {An) \ $0 > 

+ < $0 I G (Ai) X F (A2) X ....G (Ajv) I $0 > +■•• (39) 
+ < $0 I G(Ai) X G(A2) X ....F{An) I $0 > 

where F{Aj) is the operator defined as a flux surface integral associated to pixel j: 

F{Aj) = ^ [ d^s- [*1'(r')V*(r')- V*^(r')*(r')l (40) 
2mi J g L J 



'AT 



*We add the probabilities of non-exclusive events, whith provides an upper bound of the real probability of double 
detection. 
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{(Ps the differential vector perpendicular to the surface) and where G{Aj) is a volume integral associated 
to the same pixel: 

G{Aj)=[ *t(i.')^(r') (41) 

The value of J-{Ai, A2, ■■■An) provides the time derivative of the probability of detection at all selected 
pixels, which we calculate in § 12.41 

Now, because the various pixels do not overlap, the field operators commute and we can push all ^P^'s 
to the left, all ^''s to the right; then we expand these operators on a basis that has Ua and up as its two 
first vectors: 

*(r) = u„(r) a„ +u^((r) + ... (42) 

The end of the expansion, noted corresponds to the components of ^'(r) on other modes that must be 
added to modes a and (3 to form a complete orthogonal basis in the space of states of one single particle; 
it is easy to see that they give vanishing contributions to the average in state | $0 >■ The structure of 
any term in (|39|1 then becomes (for the sake of simplicity, we just write the first term): 



< 



$0 I 2^ {/a, [<(ri)4 + w;(ri)4] '^'^i • [Vw„(ri)a„ + Wup{r[)ap] - c.c.} 

N 



(43) 



where c.c. means complex conjugate and where Oa.af is the normal ordering operator that puts all the 
creation operators aj^, ^ to the left of all annihilation operators aa,p- Each term of the product inside this 
matrix element contains a product of operators that give, either zero, or always the same matrix element 
NalNpl^ For obtaining a non-zero value, two conditions are necessary: 

(i) the number of at,'s should be equal to that of a^'s 

(ii) the number of a^s, minus that of a/5's, should be equal to Na — Np. 
These conditions are fulfilled with the help of two integrals: 

d\ r dA i(N^.N^)A 

„ „ 2tt 27r ^ ' 

and by multiplying: 

(i) every UQ,(r^) by e''^, and every u*^{t'^) by e^*'^ (without changing the wave functions related to (3) 

(ii) then every Ua{v'j) by e'^, and every up{r'j) by e~*^ (without touching the complex conjugate wave 
functions). 

This provides: 



^ — f ^ gi(Af^-AfQ,)A 

5^ /^^ d^s, ■ [<(r;)Vu„(ri)e^'^ + zi^(r;)Vti^(r;)e-A + ul{v',)Vup{v\)e-^'->^+^) + ul{v\)Vu^{v',)e^'^^+^) 
-w„(r'i)V<(r'i)e»A _ up{v'^)Vu*p{v'^)e-'^ - u^{r[)S/u*p{r[)eA^+^^ - up{r[)Vul{r[)e~'(^+^'^ 



N 

Y[ [<(r;.)iia(r;)e^^ + u*p{rr)up{rr)e-'^ + ul{rr)up{r'j)e-^(^+^^ + c.c. 



J=2 



Sim. 



(45) 

where "sim." is for the iV — 1 similar terms where the gradients occur for j — 2, 3, ■■N, instead of j — 1. 

Now we assume that the experiment is properly aligned so that the wavefronts of the wave functions 
a and (3 coincide in all detection regions; then, in the gradients: 



V-ita(r) = m„(r) 



(46) 
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the vectors kc, and are parallel; actually, on each pixel j we take these two vectors as equal to the 
same constant value kA , assuming that the pixels are small and that the wavelengths of the two wave 
functions are almost equal. Then f^S]) simplifies into: 



^^r ^ r ^ e 



N„ )A h 
2m 



N 



-u^(r;-)w/3(: 



-i{\+A) , 



-iA 



-u^(r;>„(r;.)e 



J(A+A) 



(47) 

In this expression, the integrand in the A and A integrals is a product of N factors corresponding to the 
individual pixels, but this does not imply the absence of correlations (the integral of a product is not the 
product of integrals). The probability flux T contains surface integrals through the front surfaces of the 
pixels, as expected, but also volume integrals in other pixels, which is less intuitivc0. As a consequence, 
for obtaining a non-zero probability flux J- in 3A^ dimensions, it is not sufhcicnt to have a non-zero three 
dimension probability flux through one (or several) pixels; it is also necessary that some probability has 
already accumulated in the other pixels. In other words, at the very moment where the wave functions 
reach the front surface of the pixels, the first time derivative of the probability density in 3iV dimension 
remains zero, while only the A'^-th order time derivative is non-zero (this will be seen more explicitly in § 
2.3p . This is analogous to the photon detection process with N atoms in quantum optics, see for instance 
Glauber [501. 



2.3 Time dependence 

Consider an experiment where each source emits a wave packet in a finite time. We assume that, in each 
wave packet, all the particles still remain in the same quantum state, but that this state is now time 
dependent; the state of the system is then still given by but with time dependent states a and /3, 
so that the creation operators are now (oq)^ (t) and (a^)^ {t). All the calculation of the previous section 
remains valid, the main difference being that the wave functions are time dependent: Ua{r^ t) and U/5(r, t). 

We must therefore take into account possible time dependences of the wave fronts of the two wave 
functions, as well as those of their amplitudes and phases. If, for instance, the interferometer is perfectly 
symmetric, and if the two wave packets are emitted at the same time, they will reach the beam splitters 
of the detection regions at the same time with wave fronts that will perfectly overlap at the detectors; the 
amplitudes of the two wave functions will always be the same. See for instance ref. [51,52] for a discussion 
of the time evolution of the phase of Bose-Einstein condensates, including the effects of the interactions 
within the condensate. 

If we consider separately each factor inside the A and A integral of (f47|) . we come back to the usual 
three dimension space; two different kinds of integrals then occur: 



and: 

ul (r;. , t)uf,{v'^ , t)e'^(^+A) + u*p{v'^ , t)w„ (r^- , t)e'(^+A) 



(48) 



(49) 



^Equation H35| l shows that, in the definition of surface in SA'^ dimension space, the 2 dimension front surface of any pixel 
is associated with all three dimensions of any other pixel, including its depth. These dimensions play the role of transverse 
dimensions over which an integration has to be performed to obtain the flux (similarly, in 3 dimensions, the flux through a 
surface perpendicular to Oz contains an integration over the transverse directions Ox and Oy). 
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The conservation law in ordinary space implies that fj (t) is related to the time derivative of gj (t) : it gives 
the contribution of the front surface of the pixel to the time variation of the accumulated probability 
gj{t) in volume A^. The total time derivative of gj{t) is given by: 

j^9it) = fAt)-fiit) (50) 

where f^^ (t) is the flux of the three dimensional probability current through the lateral and rear surfaces 
of volume A^-; the first term in the right hand side is the entering flux, the second term the out-going 
flux, with a leak through the rear surface that begins to be non-zero as soon as the wave functions have 
crossed the entire detection volume A^. But we do not have to take into account this out-going flux 
of probability: we assume that the detection process absorbs all bosons. For instance, once they enter 
volume Aj, the atoms are ionized and the emitted electron is amplified into an cascade process, as in 
a photomultiplier; the detection probability accumulated over time does not decrease under the effect 
of fj~{t). Therefore we must ignore the second term in the rhs of ^SU\i . and replace (jH]) by the more 
appropriate definition of gj (t) : 

g,{t)= f dt' f,(t') (51) 

(we assume that time t — Q occurs just before the wave packets reach the detectors). With this relation, 
we no longer have to manipulate two independent functions / and 5; moreover, the value of gj{t) now 
depends only of the values of wave functions on the front surface of the detector, which is physically 
satisfying (while (|49p contains contributions of the wave functions in all volume A^ , an unphysical result 
if this volume has a large depth). 

We have already assumed in that the wavefronts of the two waves are parallel on every pixel; we 
moreover assume that is perpendicular to the surface of the pixel, and then call ¥'(Aj) their relative 
phase over this pixel, taking it as a constant over the pixel and over time, during the propagation of the 
wave functions (which is the case if the interferometer is symmetrical, as in the figure). Moreover, we 
assume that the two wave functions and up have the same square modulus ktAj (i) at this pixel, so 
that the interference contrast is optimal (again, this is related to a proper alignment of the interferometer). 
Then (HHI becomes: 



fj{t) ~ |kA,| |uA,(t)|^{cosA-|-cos[<^(Aj) - A - A]} = ^Pj(0 x {cos A 4- cos [(y5(Aj) - A - A]} 



with: 



(52) 

P.(i) = ^^.|kA,|^*di>A,(i')|' (53) 

where Sj is the area of pixel j. 

Finally, we assume that all the pixels are identical so that their detection areas have the same value 
5. We then obtain the simphfied expression: 

,2^ y,2^ cos[(iV^-iV„)A] -np,(t){cosA-^cos[^(A,)-A-A]} (54) 

(we have used A parity to replace the exponential in {Np — Na)A by a cosine, so that the reality of the 
expression is more obvious) and the accumulated probability at time t is: 

dX dA ^ 

— / — cos[(iV^-7V„)A] xl[pj{t) {cosA + cos [<f {A,) ^ A- X]} (55) 
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We finally consider a situation where mi pixels belong to the first detector, m2 to the second, etc., 

each sitting in one detection region after the last beam splitters. We assume that the front surface of the 
detectors are parallel to the wave fronts, so that all phases differences ip{Aj) collapse into 4 values only, 
two (in region Da) containing the phase shift (, two (in region Db) containing the phase shift 9: 

ifA — Q for the mi first measurements 

ipA — C + for the next 1112 measurements , , 

ifB + for the next 7713 measurements ^ ' 

^Pb + + -K iov the last m4 measurements 

We note that unitarity (particle conservation) requires that the third and fourth angle are obtained by 
adding tt to the first and third angles). So the time derivative of the probability of obtaining a particular 
sequence {m\,m2, ms, 1714) with given pixels is (with the new variable A' = A— A, (' = (—ipA, &' = O—fs)' 

Vit) ^ \[p,{t) j — J — cos[(7V^~7V„)A] [cosA + cos(C'-A')]™McosA-cos(C'-A')]'"= 

X [cos A + cos {0' - A')]"' [cos A - cos {0' - A')]™' 

(57) 

For short times, when the wave packets begin to reach the detectors, the probabilities Pj{t) grow 
linearly in time from zero, as usual in a 3 dimensional problem. The coincidence probability T'{t) contains 
a product of A'' values of Pj{t), so that it will initially grow much more slowly, with only a N-th order 
non-zero time derivative. For longer times, when the pj(t)'s have grown to larger values, any derivative 
of 'P{t) may be non-zero. At the end of the experiment, when the wave packets have entirely crossed the 
detectors and all the particles are absorbed, the Pj(t)'s reach their limiting value pj, and the probability 
is (from now on, we drop the primes, which just introduce a redefinition of the origin of the angles) : 

N 

P(mi,m2,m3,m4)-n^, /:^f /^^ f cos [(A^^ - 7V„)A] [cos A + cos (C - A)]"Mcos A - cos (C - A)] 

X [cos A cos (6» - A)]™' [cos A - cos {0 - A)] 

(58) 

2.4 Counting factors and probabilities 

At this point, we must take counting factors into account. There are: 

Q! 



mi!(Q — mi)! 



(59) 



different configurations of the pixels in the first detector that lead to the same number of detections mi . 
For the two detectors in Da, this number becomes: 

^ ^ (60) 



mi!(Q - mi)! m2!(Q - m2)! 
But, if we note niA = mi + m2 and use the Stirling formula, we can approximate: 
log(Q - mi)! + log(Q - m2)! 

~ (Q-mi + i) [logQ + log(l-i|^)] -(Q-mi) + (Q-m2 + i) [log Q + log (l - i^)] - (Q - m2) 

(61) 

or, if we expand the logarithms of (1 — m\^2/Q)' 

(2Q-mA + l)logQ-Q^ + ^' + ^^ -2Q + m^ + ... (62) 
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the second and the fifth term cancel each other, the third can be ignored because of ([55)): an exponentiation 
then provides the following term in the denominator of the counting factor: 



(63) 



Qm,A 

The Q! disappear, and the number of different configurations in region Da is: 

1 



g"'- (64) 
7711 !m2! 



Finally, we also have to take into account the factors in ([55)1 . These factors fluctuate among all the 
pixel configurations we have counted, since some pixels near the center of the modes are better coupled to 
the boson field and have larger 's than those that are on the sides. If we assume that the number of pixels 
of each detector is much larger than mi and 777.2, in the summation over all possible configurations of the 
pixels, we can replace each p^ by its average <p> over the detectoJl. If we assume that all detectors are 
identical, this introduces a factor < p >"^^ in the counting factor. When we take into account the other 
detection region Db, the factor Q™-*, together with the factor Q™^, can be grouped with the prefactor 
S*^ in (|57p to provide {QS)^ , which contains the total detection volume to the power iV, as naturaQ; on 
the other hand, the factor < p >^ is irrelevant, since it does not affect the relative values. At the end, 
we recover expression (fTSj) for the probability of obtaining the series of results {mi,m2, 7773, 7774). 

This calculation shows precisely what are the experimental parameters that are important to preserve 
the interesting interference effects, and expresses them in geometrical terms. The main physical idea is 
that the detection process should not give any indication, even in principle, of the source from which 
the particles have originated: on the detection surface, the two sources produce indistinguishable wave 
functions. Therefore, in practice, what is relevant is not the coherence length of the wave functions over 
the entire detection regions, as the calculation of §[T] could suggest, since in these regions the modes are 
defined mathematically in a half-infinite space; what really matters is the parallelism of the wave fronts 
of the two wave functions with the input surface of the detectors. Moreover, if necessary, formulas such as 
(011) and (|i5| allow us to calculate the corrections introduced by wave front mismatch, and therefore to 
have a more realistic idea of the experimental requirements; for instance, if the and are not strictly 
parallel and perpendicular to the surface of the photodetectors, one can write write ha.p = ± <5k(r) 
and calculate the correction to (05]) to first order in 6k, etc. 



3 EPR argument and Bell theorem for parity 

EPR variables are pairs of variables for which the result of a measurement made by Alice can be used to 
predict the result of a measurement made by Bob with certainty. For instance, the numbers of particles 
detected by Alice and by Bob are such a pair, provided we assume that the experiment has 100% efficiency 
(no particle is missed by the detectors): Alice knows that, if she has measured niA particles. Bob will 
detect — 777^ particles. It is therefore possible to use the EPR argument to show that niA and 777^ 
correspond to elements of reality that were determined before any measurement took place. Moreover, 
this also allows us to define an ensemble of events for which ttia and ttib are fixed as an ensemble that 
is independent of the settings used by Alice and Bob; this independence is essential for the derivation of 

''If mi = 1, the summation provides exactly < p > the average by definition. If mi = 2, the second pixel can not coincide 
with the first, so that the average of the product pip2 is not exactly < p >^; nevertheless, if the number of pixels Q is 
much larger than 2, the average is indeed < p >^ to a very good approximation. By recurrence, as long as the number of 
detections m remains much smaller than the number of pixels, one can safely replace the average of the product by the 
product of averages. 

'^The probability remains invariant if, at constant detection area, the value of the number of pixels Q is increased. 
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the Bell inequalities within local realism [10]. So we may either study situations where tua and ttib are 
left to fluctuate freely, or where they are fixed (as with spin condensates [16]). 

When Na ~ Np, we have seen in ? ll.3.2l that another pair of EPR variables is provided by the parities 
A = (—1)™^ and B = (—1)™* of the results observed by Ahce and Bob: if they choose opposite values 
( = —9 for their settings, perfect correlations occur, even if Alice and Bob are at an arbitrarily large 
distance from each other. We now study quantum violations of local realism with these variables. 



3.1 Parity and BCHSH inequalities 

We suppose that in the experiment of Fig. 1, Alice and Bob each use two different angle settings, ( 
and for Alice and 9 and 9' for Bob. Within local realism (EPR argument), for each realization of the 
experiment the observed results depend only on the local settings. We can then define A as the parity 
observed by Alice if she chooses setting C,, and A' the parity if she chooses similarly. Bob obtains 
results B or B' depending on his choice 9 or 9'; all these results are parities equal to ±1. Then, since 
either B + B' oi B — B' vanishes, within local realism we have the relation: 



2>AB + AB'+A'B - A'B' < 2 



(65) 



For an ensemble of events, the average of this quantity over many realizations must then also have a value 
between -2 and +2 (BCHSH theorem). 

In quantum mechanics "unperformed experiments have no results" [53] : any given realization of the 
experiment necessarily corresponds to one single whole experimental arrangement, and it is never possible 
to define simultaneously all 4 numbers in Eq. (j65p . One can nevertheless calculate the quantum average 
of the product of the results for given settings, and derive the expression: 



Q = (AB) + {AB') + (A'B) - (A'B'] 



(66) 



but there is no reason to expect Q to be between —2 and +2. 

Since ([24|) reduces to (|15p (with M = N) when R = and T = 1 , we can proceed from the more 
general formula ([M]) . The calculation of the average (AB) is very similar to that of section (iv) of the 
Appendix, but with a factor (— l)™2+m4 included in the sum on mi, - • • ,7714. The equivalent of (|107p . 
obtained after summations over mi and m2 (with constant sum m^) and over and nii (with constant 
sum ms) is: 



M T}N-M r+'^ dA 



/_7fAcos[(A^„-iV^)A] [cos A] 
/:^|A[2cos(C + A)]'"- [2cos{9-X)Y 



N-M 



(67) 



Formula p03p of the Appendix can then be used, with M replaced by iV — M . Therefore we see that the 
average {AB) of the product vanishes, unless the two following conditions are met: 



M is even 

M < 2Na and M < 2Np 
It these two conditions are met, the first line of ((H7|) becomes: 



m^lms! (7V„-f )!(iV^-f)! 
while the second line provides, with the help of formula of the Appendix: 



M 



dX 

2^ 



cos A + 



c- 



M 



Ml 



M 



M 



(68) 



(69) 



(70) 
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Finally, the sum over tua and ms (with constant sum M) gives the result: 



NJNf3\ 



_rpM J^N~M 



COS 



M 



(71) 



If M is left to fluctuate, a summation of this expression over M should be done. But another point 
of view is to decide to count only the events where M is fixe41- Then this average should be compared 
with the probability that M particles will be detected, given by formula (|109p in the Appendix; dividing 
((7T|) by (fTM)) now provides: 



{AB) 



NJNplMliN - M)\ 



m (iV.-f)!(iV,-f)![(f)! 
In the case M = N, the second condition ([68| requires than Na 



M 



(AB) = cos 



N 



for M = N 



(72) 

Np — M/2, in which case we get: 

(73) 



One can put this into Q of Eq. ([55]) : Alice's measurement angle is taken for convenience as (j)a = 2^ 
and Bob's as — —29. Then defining E{(j)a — 4>b) — cos^ {<j)a — 4>b) and setting (jj^ — 4>b — (pb — 4>a' = 
0f)' ~ (pa — ^ and (pb' — 4'a' ~ Slu we can maximize Q — 3E{uj) — E(Suj) to find the greatest violation 
of the inequality for each N. For = 2 we find Qmax = 2.41 in agreement with Ref. [37]; for A^ = 4, 
Qmax — 2.36; and for A^ — > oo, Qmax ^ 2.33. These values are obtained for a value of the angles 
corresponding to ui = y/ln3/N, which decreases relatively slowly with A^. The conclusion is that the 
system continues to violate local realism for arbitrarily large condensates. As already noted in § 11.3.11 
this is a direct consequence of the effects of the quantum angle A, since no such violation could occur if 
this angle was zero. 

Suppose now we measure M = N — 1 particles with A^^ = M/2, Na = M/2 + 1. Then the coefficient 
of the cosine in Eq. ([72]) is (M/2 + 1)/(M + 1), which is 2/3 at M = 2 and smaller for larger M, so that 
this case never violates the BCHSH inequality since 2/3 x 2.41 — 1.61 < 2. If ([7T|) had been used instead 
of l|72p . we would be even even further from any violation, since the first average value is smaller than 
the second. The conclusion is that even one single missed particle ruins the quantum violation. 



3.2 Three Fock states and three interferometers; GHZ contradictions 

With a triple- Fock state source (TFS) as shown in FigHlwe can demonstrate GHZ contradictions [13, 14]. 
Such a contradiction occurs when local realism predicts a quantity to be, say, +1 while quantum mechanics 
predicts the opposite, —1. Previous such contradictions were carried out with states known variously 
as GHZ states, NOON states, or maximally entangled states. These wave functions are of the form 

u|-|--|--t----)+i;| •••) with particular values of the phases u and v. The original GHZ calculations 

[13] were done with three- and four-body NOON states, and this was generalized to A^ particles by 
Mermin [12]. Yurke and Stoler [38] showed how an interferometer with three one-particle sources also 
could give a GHZ contradiction. We will replace their sources with Bose condensates to show how new 
Af-body contradictions can be developed. 
The initial TFS is: 

1$) ^ \Na,N,,N,) ^ ^^^,^^,^^, 4^°«r«r |0) (74) 

®With the experimental setup of fig. [2] one can include the results of measurements given by the detectors in channels 
5 and 6 in the preparation procedure; only the events in which ms -|- me = N — M are retained in the sample considered. 
Since this procedure remains independent of the settings 6 and this does not open a "sample bias loophole". 



19 




Figure 4: Interferometer with three Fock-state condensate sources and three detectors. The particles 
from each source can reach two detectors. Each detector has two subdetectors, which wiU register a +1 
for the odd- numbered subdetector and —1 for the even- numbered. We average the quantity ABC where 
A = ±1 for Ahce's detector, and B — ±1 for Bob's, and C = ±1 for Carole's. 



As in § [T] the output modes (destruction operators) ai ■ ■ ■ can be written in terms of the modes at the 
sources a^, and with three phase shifts of ^, x, or 0. We find: 



ai = - [e'"^aa - iap] 



a2 



03 

as 



lie"' an 



"4=2 [*e**a/3 - a^] , 



(75) 



We write generally ai = ViaUa + Vipap + Vi^a^. We consider only the case where every particle in the 
source is detected, so the probability that we find rrii particles in detector i = 1 • • • 6, is: 

1 



TOl! • • • JTlg! 



|(0|ar...ari^a,^/3,^7>r 



(76) 



(this relation is actually an equality, but we write it only as a proportionality relation since we will change 
the normalization below). We can develop the matrix element just as we did in § 11.11 



i=l 

mi me 

E-E 

Pl = P6 = 



n 



Pia, Pifi Pi-I 



(77) 



where the sums are over all pia, pip and pi^ such that pia +Pif3+Pi'y = rrii. We now replace the (5- functions 
by integrals: 

with similar integrals over Xp and A-y. In the sum above then, we have every replaced by (fiae'^°)^"* 
etc. so that we can redo the sums over the pic, etc. to find the probability for the rrii arrangement under 
the condition that all the source particles are detected: 

P(mi, • • • , me) : / dr' f (ire-*[^=(^°-^=)+^^(^^-^;^)+^-(^--^^)l 17 {n'*n,r' (79) 



20 



where $7^ = Viae^^" + Viije^^i^ + Vi^e^^i and Vl'^ has the same expression with primed A's; dr represents 
the integrals over A^, Xp, and A-y, and dr' over the A^, A^, and A^^. 

In a ideal experiment with 100% detection efficiency, the numbers of particles detected in each region 
are EPR variables, since the value of two of these variables determines the value of the third with certainty; 
these perfect correlations are independent of the settings (phase shifts of C, X, or 6'), so that choosing the 
number of detections in each region defines a class of events that is independent of the settings. Here, 
assuming that each source emits 7V/3 particles (otherwise, we find zero average values, see below): 

we will also assume that each detector registers exactly iV/3 particle^. We can put in this restriction, 
when we sum on mi • • • mg to get averages, by including three (5-functions of the form: 



2n 



(80) 



with similar ones specifying + — N/S and + mg = N/3. The rrii sums are then done 
independently of one another giving a normalization sum of: 



Af= dr, dr' / dre-'^[^=(^=-^»)+^^(^^-^^)+^.(^.->;)l 



^ ^-^N/^pA+PB+Pc] exp 

where = p2 ^ pA, Ps ^ Pi ^ Pb, and p^ ^ pe ^ 
The sum in the exponential is easily done: 

1 



,tpi 



(81) 



PC and J dTp represents the new three-fold integration. 



1=1 

^g-«(^/3-A^) (gipB _|_ gJpA) _|_ gj(A-,-Al,)] (gipc _|_ gVB) 



(82) 



We expand the exponential of this quantity in series in e e and e^^^'' ^t-* and do the 

integrals. When each source emits exactly particles, we obtain: 



1 



1=0 



(83) 



Similarly we average the quantities A, B, and C each ±1 measured by Alice, Bob, and Carole according 
to: 

{ABC)= J2 ,mg) (84) 

mi ■ ■■rriQ 

where the prime on the sum means we again restrict the sums to the case of iV/3 particles reaching each 
detector. With this requirement the average vanishes unless each source emits exactly iV/3 particles, 
which is why above we considered just that case. We have then: 



jV/3! 



(ABC) = 



iN/3-qy.q{ 



jV/3! 



q V(W/3-g)!g! 



(85) 



^We have also performed more general calculations — not given here — in which this restriction does not apply, but then 
we have found no GHZ contradictions. 
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Note that if(' + 6' + x = 0wc find (ABC) = 1: perfect correlations exist between the resuhs since their 
product is fixed. Thus, if we know the parity of the results of two of the experimenters, we immediately 
know that of the third, even if that person is very far away. Thus an EPR argument applies to these 
variables. 

For the case = 3, we find (ABC) = cos{( + d + x), which is the same form of the original GHZ case 
found from a three-body NOON state. This result agrees with the interferometer result of Ref. [38] as 
expected. Local realism predicts that, for each realization of the experiment, the product of the results 
is given by a product A(QB(9)C{x)- To get agreement with quantum mechanics in situations of perfect 
correlations we must have: 

A(7r/2) B{tt/2) C(0) = -1 
A(7r/2) B{0) C{tt/2) ^ -1 

A{Q) B{tt/2) C{tt/2) = -1 (86) 

But then we obtain by product A{Q)B{0)C{0) = —1, while quantum mechanics gives +1, in complete 
contradiction. In our case we get new contradictions for larger N; consider for instance = 9, in which 
case: 

{ABC) = ^ [27 cos(C + + x) + cos 3(C + ^? + x)] (87) 

The above argument goes through exactly in the same way. More generally, any time N/ 3 is odd we get 
a similar result for arbitrary TV. 

Thus the TFS provides new GHZ-type contradictions for TV particles without having to prepare NOON 
states. 

3.3 Hardy impossibilities 

Hardy impossibilities are treated by use of the interferometer shown in Fig. [Sj based on the one dis- 
cussed in Ref. [46] for N ^ 2. The heart of the system is the beam splitter at the center; due to Bose 




Figure 5: An interferometer with particle sources a and f3, with beam splitters designated by BS and 
mirrors by M. In both detection regions, the detectors at may be replaced by the D^, placed before 
the beam splitters. 
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interference it has the property that, if an equal number of particles approaches each side, then an even 
number must emerge from each side. The detection beam splitters BSA and BSB are each set to have 
a transmission probability of 1/3, and the path differences are such that, by destructive interference, no 
particle reaches if only source is used; similarly, no particle reaches D3 if alone is used. Alice 
can use either the detectors D12 after her beam splitter, or D'^ 2 before; Bob can choose either 03^4, or 
D3 4. This gives 4 arrangements of experiments: DD^ DD' , D'D, or D'D', with probability amplitudes 
Cxy (tii, m2; m3, 7714), where XY is any of these 4 arrangements and the m values are the numbers of 
particles detected at each counter. 

We find the destruction operators for the detector modes as we have done in previous sections. For 
the primed detectors we find: 

and for the unprimed detectors: 

aD^=-^a^ + ^a0, an, = -^a^ ^^^^ 

In general we write these results as: 

di = ViaO-a + Vipap (90) 

Note that, because of the 1/3 transmission probability at BSA and BSB, Bose interference causes Oq and 
ap to drop out of the second and third of Eqs ([SQ]) . respectively. 
The amplitude is given by: 

6 

C'xY{mi,m2;m3,m4) ^ (0| (uiaaa + Wj/ja/j)™' ajj^^aj,^" |0) (91) 

2=1 

As we have done in previous sections, we expand the binomials, evaluate the operator matrix element, 
replace the resulting (5-functions by integrals and resum the series to find: 

CxY(™i,m2;m3,m4) - T ^ T ^e"'""^" e"'"^^^ n i-^^e^'^ + v^pe^'^p (92) 

In all the following wc assume that each source emits A^/2 particles, where N/2 is odd, and that 
detector A and detector B each receive exactly A^/2 particles; this is possible since, as above, the number 
of particles detected in each region can define a sample of realizations that is independent of the settings 
(we have to make this assumption since it turns out that the argument works only in this case). First 
consider both Alice and Bob using primed detectors. The amplitude for receiving A^/2 particles in each 
of D2 and Dg is: 



This quantity must vanish because N/2 is odd. This situation is an example of the beam splitter rule 
mentioned above. The result is that and D3 cannot collect all the particles if N/2 are detected on 
each side. 
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Consider next the case where one experimenter uses a primed set of detectors and the other the 
unprimed: 

~ Sm'.fi (94) 

The factor (e*'^'' ) combines with the first exponential so that 7714 must vanish. This quantity vanishes 
because of the destructive interference effects at BSA and BSB caused by the 1/3 transmission probabihty 
of the beam splitters; but Cdd'(Oi t! T'^) ^- Thus, if Alice observes N/2 particles at D2, when Bob 
uses the primed detectors he observes with certainty N/2 particles at D3; similarly, if Bob has seen N/2 
particles in D3, in the D'D configuration Alice must see N/2 in . 

Wc now consider events where both experimenters do unprimed experiments and each of them finds 
N/2 particles in D2 and D3; the corresponding probability is: 

Coo(0,f ;f ,0)~£^ £^^-^f(^..V (,...)-/^(,..„)-/Vo, (95) 

(for N = 6, the normalized value is 1/216), which means that events exist where A''/2 particles are 
detected at both detectors D2 and D3. However, in any of these events, if Bob had at the last instant 
changed to the primed detectors, he would surely have obtained A^/2 particles in D3, because of the 
certainty mentioned above (while Alice still has N/2 particles in D2). Similarly, if it is Alice who chooses 
primed detectors at the last moment, she always obtains N/2 particles in (while Bob continues to 
have N/2 particles in D3). Now, had both changed their minds after the emission and chosen the primed 
arrangement, local realism implies that they would have found A''/2 particles each in D2 and D3: such 
events must exist. But the corresponding quantum probability is zero, in complete contradiction. The 
result is the Hardy impossibility of Ref. [46] generalized to particles. 



4 Conclusion 

Fock-state condensates appear as remarkably versatile, able to create violations that usually require 
elaborate entangled wave functions, and produce new AT-body violations. Compared to GHZ states or 
other elaborate quantum states, they have the advantage of being accessible through the phenomenon of 
Bose-Einstein condensation, with no limitation in principle concerning the number of particles involved. 
By contrast, the production of GHZ states requires elaborate measurement procedures, so that it seems 
difficult to produce them with more than a few particles (to our knowledge, the present world record 
is 5, see [54]); moreover, they are much more sensitive to decoherence, which destroys their quantum 
coherence properties [33]. 

From an experimental point of view, the major requirement is that all particles present in the initial 
double Fock state should be detected, which will of course put a practical limit on the number of particles 
involved. Using Bose condensed gases of metastable He atoms seems to be an attractive possibility, since 
the detection of individual atoms is possible with micro-channel plates [18,19]. With alkali atoms, one 
could also measure the position of the particles at the outputs of interferometers by laser fluorescence, 
obtaining a non-destructive quantum measurement of mi, ..1714. The realization of interferometers seems 
also possible, since intcrferomctry with Bose-Einstein condensates has already been performed [55] with 
the help of Bragg scattering optical beam splitters [56,57]. Another possibility may be to use cavity 
quantum electrodynamics and quantum non-demolition photon counting methods [58] to prepare multiple 
Fock states. Experiments therefore do not seem to be out of reach. 

Laboratoire Kastler Brossel is "UMR 8552 du CNRS, de TENS, et de I'Universite Pierre et Marie 
Curie". 
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APPENDIX I 

In this appendix, we give some formulas that are useful for the calculations of this article, in particular 
to check the normalization of (fT5|) and (|24|) . 

(i) Wallis integral. We consider the integral: 

K = [cosXf (96) 

which, if the limits are changed to and 7r/2, becomes a Wallis integral (divided by 2tt if one takes the 
usual definition of these integrals). Expanding the integrand provides: 

2-^ [e^^ + e-^^l = 2-^ y „ e^^^'^^)^ (97) 

The only exponentials that survive the A integration are those with vanishing exponent (N — 2q = 0). 
Therefore: 

if N is odd, K ^0 

if TV is even, if = 2-^ (98) 

(ii) Normalization integral. We define: 

J = ^ cos [{N^ - Np) A] [cos A]'^ = Re ( f^^ ^e^(^»-^^)A [cos A]'^| (99) 

and expand: 



[cosAf 



-, M M 



2 

Only a term with M ^ 2q — Na — Np can survive the integration, so that: 



^g!(M-g)! 



(101) 



2 

Therefore, J is non-zero only if: 

M has the same parity a.s N^- Np ; - M < Na~ Np < +M (102) 

and then: 

J = 2-^-. ^ (103) 

(iii) Normalization of We now consider the probabilities V {mi, 1712,1^3,1714) given by ([T5)) and 
calculate their sum over mi, m2, TO3, 7714, when these variables have a constant sum N. We do these sums 
in three steps: a summation over mi and TO2 (with constant sum uia), a summation over 1713 and 7714 
(with constant sum m^), and a summation over uia and ms (with constant sum M). The first two sums 
reconstruct powers of a binomial, the A integral disappears, and we obtain: 

N^lNpl 2-^ ^ cos [{N^ ~ Np) A] x [2cosA]^ (104) 

which, with pjl^jl for M ^ N = Na+ Np, gives: 

2-^ X (105) 
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Then the summation over tua and ms gives: 

2-^(1 + 1)^ = 1 (106) 

as expected. 

(iv) Normahzation of ([24| . We now consider the probabilities P(toi, m2, TO3, 7714) given by ([24| and 
calculate their sum over any toi, m2, 7713, 7714. We do this by the same three summation as above (with 
77^^ +771B = M, instead of N), plus a summation over M ranging from to iV. The first two summations 
give: 

AT I AT I nN — 2M r+TT j\ 

f UN A/n. ^"'^"'''' / i^cos[(iV„-7V,)A][cosA]^-^-^[2cosA]^-^ (107) 
mAmBi (N - My. J_^ 2tt 

or, when (|103p is inserted: 

ATI n-Af 



77iA!mB! (iV- M)! 
The summation over ttia and 777,3 with constant sum AI then gives: 

iV! 



M! {N - My. 



which provides the probability of detecting M particles, independently of which of the 4 detectors is 
activated. This probability is maximal when: 

N-MT M , , 

-M-R-' or -^T (110) 

The larger the transmission coefficient T, the larger the most likely value of M, as one could expect 
physically. Finally, a summation of (|109p over M between and A'' gives: 



{R + Tf = 1 (111) 



and the total probability is 1, as expected. 



APPENDIX II 

In this appendix, we investigate how Eq. (|15p is changed when the initial state |$o) is different from 
the double Fock state considered in ([1]). 
(a) Coherent states 

We first assume that each of the modes a, /3 is in a coherent state with phases 4>a, 4>i3 and the same 
amplitude E: 

|$o) = ® 10/3) (112) 
with the usual expression of the coherent states: 

J^\Nc.,p=r) (113) 

The calculation of § 11.11 is then simplified since this state is a common eigenvector of both annihilation 
operators Ua and a^. There is no need to introduce conservation rules, and neither A nor A enter the 
expressions. Eq. p5|) becomes: 

V {mi, 1712, ms, 1124) [1 + cos(C + (f>a ~ 'f'p)]"'^ [1 - cos (C + (t)a- 0/3)]"^ /-, -, .X 

X [1 + cos (e + ^fj- ^^)r' [1 - cos i-e + (bo^T' ^ ' 
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Here, no A distribution occurs, in contrast with the relative phase of the two states is perfectly 
defined and takes the exact value (j)^ — (fif}. Now, we can also assume that the initial phases of the two 
coherent states completely random. Then, an average over all possible values of 0ctand (f>/3 leads to: 

r{rm,m,, mg, m^) ^ J ^ [1 + cos (C + 0)]'"^ [1 - cos (C + 0)]"^ , 

[1 + cos i-e + 0)]"^ [1 - cos {-0 + 0)]™^ ^ ' 



We now obtain an expression that is similar to Eq. ()15p . but with a difference: the terms of the product 
in the integral are always positive, as if the quantum angle A had been set equal to zero; no violation of 
Bell inequalities is therefore possible. This was expected: with coherent states, the phase pre-exists the 
measurement and is not created under the effect of quantum measurement, as was the case with Fock 
states; an unknown classical variable does not lead to violations local realism. Moreover, the requirement 
of measuring all particles does not apply in this case, since the initial state does not have an upper bound 
for the populations. 

(b) Phase state 

We now assume choose a state that has a fixed number of particles, but a well defined phase (f>o 
between the two modes: 

\<^o,N)^-l= [e^^o,i+al]"\0)^VN^.Y: -^7^ (4)" (4)"" |0) (116) 



with N even. We assume that all particles are measured: ^ ■ rm = N. The probability amplitude we 
wish to calculate is: 

1 4 

C™,...™, = (0| n a7^- |$o, iV) (117) 

3 = 1 

where the are defined in Q and written more generically in ([!]). We will use two different methods 
to do the calculation, first a method based on the specific properties of phase states, and then a more 
generic method extending the results of § 11.11 

(i) The phase state |$0i-^) is by definition a state where all bosons are created in one single state 
[e'"^" |a) + 1/3)] /\/2, none in the orthogonal state [-e*"^" |a) + /\/2. Therefore the action of the two 
annihilation operators: 

(118) 

is straightforward: the former transforms |$o, A^) into |$07 ^ ^ 1)7 the latter gives zero. Now, we can use 
([U and (jllSp to expand each ai as: 

a/2 a/2 
where the action of a^^^^ gives zero. We conclude that: 

aj \^o,N) = \[^{e''^''vjc. + v,p) |$o, N - I) (120) 

in which case: 
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The probability is then: 

1 1 7 3 ' A — 1 

W (122) 

X [(1 + cos(0 - 0o)]"-M(l - cos(0 - 0o)]"^ 

As in case (a), we have a state for which the quantum angle vanishes so that no violation of the BCHSH 
inequalities can take place. 

(ii) We can also do the calculation by a method that is similar to that of 11.11 From (|116p and (|117p . 
we obtain: 

The calculation is the same as in § II. 1[ with replaced by q and TVg by — g; the prefactors of ((T]) 
and (|116p combine to introduce a factor N\ / q\ [N — q)\, and the equivalent of pU)) is now: 



/ l^e''"-^-'" n (-..e- + ...e--)™' (123) 



In the probability, a sum over q and g' appears, including term q ^ q' corresponding to non-diagonal 
probabihty terms between two different states {N^ — q' ; Np — N — q'\ and \Na — q ; Np — N — q). One 
finally obtains: 

7^(mi,m2,m3,m4)^ , f , 2'^ T ^ T ^ G (A, A) [cos A + cos (C + A)r 

mi\m2\m3\mi\ J_^2ttJ_^2ii (124) 

X [cos A - cos (C + A)]"^ [cos A + cos {9 - A)]'"=' [cos A - cos {9 - A)]"" 

with: 

TV 

G(A,A)= gi^-l'Y't'o-X) ^r{N~q~q')K ^25) 

q,q'=0 

The result is therefore similar to (fT^ . except for the presence of the function G(A, A), which introduces 
a distribution of the phase A and of the quantum angle A. 

Equations p24p and p25p are equivalent to (|122p , although they do not contain the same distribution 
G(A,A). Equation p22p corresponds to an infinitely narrow distribution, since it can be obtained by 
replacing in (|124p G (A, A) by the product 5{X — 0o) 5{K); by contrast, (|125p defines a distribution with 
finite width. This illustrates the fact that, in (|124p . different distributions G (A, A) may lead to the same 
set of probabilities. 

(c) General state 

Consider finally the more general state |<i>o) combining two modes with a fixed total number of particles 
can be written as: 

N 

\<^o)=Y,x,\N^^q-Np = N-q) (126) 

where the complex coefficients Xq are arbitrary. The calculation is similar to that of §(ii) above, but now 
one obtains p24p with a different expression of G(A, A): 

G(A, A) = ^ XqX^, e*(^-9-9')e-«A(g-g') ^qn,^N -q')\q\{N -q)\ (127) 
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Depending on the choice of the coefficients Xq, one can build states in which the initial phase is well 
determined, as in § (ii), or completely indetermined as for Fock states; a similar conclusion holds for the 
quantum angle A. 
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